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Universal behavior of crossover scaling functions for continuous phase transitions 
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We consider two different systems exhibiting a continuous phase transition into an absorbing 
state. Both models belong to the same universality class, i.e., they are characterized by the same 
scaling functions and the same critical exponents. Varying the range of interactions we examine 
the crossover from the mean-field-like to the non-mean-field scaling behavior. A phenomenological 
scaling form is applied in order to describe the full crossover region which spans several decades. 
Our results strongly supports the hypotheses that the crossover function is universal. 
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The critical behavior of a system exhibiting a second 
order phase transition with non-mean-like scaling behav- 
ior is strongly affected by the range of interactions. The 
longer the range of interactions the stronger will be the 
critical fluctuations reduced. In the limit of infinite inter- 
actions the system is characterized by the mean-field scal- 
ing behavior. But according to the well known Ginzburg 
criterion fl|, mean-field- like behavior occurs even for fi- 
nite interaction ranges sufficiently far away from the crit- 
ical point. A crossover to the non- mean- field scaling be- 
havior takes place if one approaches the transition point. 
Although crossover phenomena are well understand in 
terms of competing fix points of the corresponding renor- 
malization group approaches (see for instance 2]), some 
aspects of crossover phenomena are still open. For in- 
stance it is an open question whether the so-called ef- 
fective exponents fulfill certain scaling relations over the 
entire crossover region (see 0, S IM^Q and references 
therein). A second open question is the main theme of 
this paper and addresses the universality of the crossover 
scaling functions. The range where the universal critical 
scaling behavior applies is usually restricted to a small 
vicinity around the critical point. Therefore it is ques- 
tioned that the full crossover region, that spans several 
decades in temperature or conjugated field, can be de- 
scribed in terms of universal scaling functions. Renor- 
malization group approaches predicted a non-universal 
behavior if one uses finite cutoff lengths whereas infinite 
cutoff lengths (which corresponds to an unphysical van- 
ishing molecular sizeMead to a universal crossover behav- 
ior (see for instance |2lla|)- On the other hand the exper- 
imental situation is also unclear since measurements over 
the whole crossover region are difficult and accurate re- 
sults are rare (see '3| for a short discussion) . Thus several 
attempts were performed in order to address this ques- 
tion via numerical simulations. For instance the two- and 
three-dimensional Ising model with various interaction 
ranges is considered in a series of papers Us- 
ing a sophisticated cluster algorithm for long-range inter- 
actions it was possible to cover the full crossover region. 
In particular a collapse of the susceptibility for different 



values of the interaction range was observed. Thus the 
crossover can be described by a single scaling function in 
agreement with renormalization group approaches. But 
this result docs not present an evidence that this scaling 
function is universal, since only one system of a given 
universality class was considered. 

The purpose of this paper is to demonstrate via numer- 
ical simulations that the crossover from non-mean-field 
to mean-field-like scaling behavior can be represented by 
universal functions. We therefore consider two different 
systems exhibiting a continuous phase transition; both 
belong to the same universality class. The dynamics of 
the models is characterized by simple particle hopping 
processes, i.e., various interaction ranges can be easily 
implemented and highly accurate data are available. In 
this way it is possible to observe the full crossover region. 
Notice that we focus in our investigations on the particu- 
lar universality class of absorbing phase transitions only 
for technical reasons. The demonstrated universality of 
crossover scaling functions can be applied to continuous 
phase transitions in general. 

The first considered model is the so-called conserved 
lattice gas (CLG) which was introduced in In the 

CLG lattice sites may be empty or occupied by one par- 
ticle. In order to mimic a repulsive interaction a given 
particle is considered as active if at least one of its neigh- 
boring sites on the lattice is occupied by another particle. 
If all neighboring sites are empty the particle remains in- 
active. Active particles are moved in the next update step 
to one of their empty nearest neighbor sites, selected at 
random. 

The second model is the so-called conserved transfer 
threshold process (CTTP) [ni. Here, lattice sites may 
be empty, occupied by one particle, or occupied by two 
particles. Empty and single occupied sites are considered 
as inactive whereas double occupied lattice sites are con- 
sidered as active. In the latter case one tries to transfer 
both particles of a given active site to randomly chosen 
empty or single occupied nearest neighbor sites. 

In our simulations (see for details) we have used 

square lattices of linear size L < 2048. Every simulation 
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starts from a random distribution of particles. After a 
transient regime both models reach a steady state char- 
acterized by the density of active sites The density p„ 
is the order parameter and the particle density p is the 
control parameter of the absorbing phase transition, i.e., 
the order parameter vanishes at the critical density Pc 
according to, p^ cx Sp^ , with the reduced control param- 
eter Sp — p/ Pc — I- Additionally to the order parameter 
we consider its fluctuations Ap^. Approaching the tran- 
sition point from above (Sp > 0) the fluctuations diverge 
according to (see IHlia) Ap, (X dp-'''. Below the criti- 
cal density (in the absorbing phase) the order parameter 
as well as its fluctuations are zero in the steady state. 

It was shown recently that the order parameter as well 
as its fluctuations obey the scaling forms Q 

p^{6p,h) ^ R{a^Sp\,a,h\''), (1) 

a^Ap^{Sp,h) - X''' D{a^SpX,a^hX''), (2) 

where h denotes an external field which is conjugated to 
the order parameter |0| . The universal scaling functions 
i?(x, y) and D{x, y) are the same for all systems belong- 
ing to a given universality class whereas all non-universal 
system-dependent features (e.g. the lattice structure, the 
update scheme, etc.) are contained in the so-called non- 
universal metric factors a^, a^, and Oa jlSj. The uni- 
versal scaling functions are normed by the conditions 
^(1,0) = -R(0,1) = i)(0,l) = 1 and the non-universal 
metric factors can be determined from the amplitudes of 
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These equations are obtained by choosing in the scaling 
forms [Eqs. 10 13)] a^SpX — 1 and a^/iA'^ = 1, respec- 
tively. 

Usually scaling functions are only known above the 
upper critical dimension where the mean-field the- 
ory applies. In the case of the CLG model and CTTP 
the mean-field scaling functions are given by 0, 
RuF{x,y) = x/2 + {y+{x/2)'^y/'^ as well as b„p{x,y) = 
R-Mp{x,y)/{y + (a;/2)^)^/^ i.e., the mean-field exponents 
are Pm.f = 1, CTm^- = 2, and = (corresponding to a 
finite jump of the fluctuations). Below the upper critical 
dimension the universal scaling functions depend on the 
dimension and are unknown due to a lack of analytical 
solutions. 

In the original CLG model and the original CTTP par- 
ticles of active sites are moved to nearest neighbors only, 
i.e., the range of interactions is i? = 1. In the following 
we consider a modified CLG model and a modified CTTP 
where particles of active sites are moved (according to the 
rules of each model) to randomly selected sites within a 
radius R. The order parameter is plotted in Fig.^for var- 
ious ranges of interactions {R G {1,2,4,..., 128}). In the 
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FIG. 1: The order parameter of the CTTP for various 
values of the interaction range R. With increasing interac- 
tion range the critical density tends to the mean-field value 
PcR-'oa = 1/2. The inset displays the order parameter which 
is rescaled according to Eq. (1181 . The dotted line correspond 
to the scaling laws y = mx with m = (ap__R=i/Pc,H=i)''°=^. 

following we examine how the varying interaction range 
affects the scaling behavior in the vicinity of the absorb- 
ing phase transition which now takes place at the critical 
density 

The crossover scaling function at zero field has to incor- 
porate the range of interactions as an additional scaling 
field. We make the phenomenological ansatz 

p.(p,i?„„) ^ X-^^- m{a,{p- p^.^) X,a-'R:J X^), (6) 

where the scaling function 9\ is universal since we al- 
low for the non- universal metric factors Up and Or. The 
Ginzburg criterion states that the mean-field picture is 
self-consistent in the active phase as long as the fluc- 
tuations within a correlation volume are small com- 
pared to the order parameter itself (see 0|). Thus the 
crossover exponent is given by = (2/?Mi? — v^fD)/ D = 
(4 — D)/2D, where i^mf = 1/2 denotes the critical expo- 
nents of the spatial correlation length. In order to avoid 
lattice effects we use the effective interaction range Q 

= lY.\^.~rf, \r^~r^\<R (7) 

where z denotes the number of lattice sites within a ra- 
dius R (see Table HJ. The mean-field scaling behavior 
should be recovered for R ^ oo, thus 

m{x,0) = 7?mf(2:,0) = a;'^"^ (8) 

which implies Up = ap^R^^/ p^^R^^ These factors were 
already determined in previous works where absorbing 
phase transitions with infinite particle hopping were in- 
vestigated jl6|. The non- universal metric factor an has 
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TABLE L The range of interactions R, the corresponding 
number of next neighbors z on a square lattice and the effec- 
tive range of interactions R^it for which we have carried out 
simulations. Additionally, the values of the critical densities 
are listed. 



R 



Zo = 2 Rcif,n = 2 



CLG 



1 

2 

4 

8 

16 

32 

64 

128 



12 

48 

196 

796 

3208 

12852 

51432 



1546 

49 
25274 

199 
204875 

401 
13146247 

6426 
105255421 



0.34494(3) 

0.22432(4) 

0.16802(7) 

0.14050(9) 

0.12977(10) 

0.12598(11) 

0.12499(16) 

0.12465(19) 



0.69392(1) 
0.63649(2) 
0.55005(3) 
0.51688(4) 
0.50552(6) 
0.50161(7) 
0.50046(8) 
0.50019(9) 



to be determined by a second condition. Several ways are 
possible (e.g. $K(0, 1) = 1) but for the sake of convenience 
we force that 9\ scales as 



for 



0, 



(9) 



where denotes the non-mean-field order parameter ex- 
ponent of the corresponding Z?-dimensional system. Set- 
ting a~^R~:i} A*^ = 1 in Eq. © yields for zero field 

(10) 

Taking into account that the Z?-dimensional scaling be- 
havior is recovered for i? = 1 we find 



^p,R=l Pc,R- 



4>I3o/{0mf-I3d) 



(11) 



According to the above scaling form we plot in Fig.[21 
the rescaled order parameter pa (aai^eff)^ as a function 
of the rescaled control parameter Op(/9 — Pc,fl)(aR-Rotf)^ 
for the two-dimensional {(p = 1/2) CLG model and the 
two-dimensional CTTP. The values of the metric factors 
are listed in TablellTl and are determined from data of 
previous simulations (via a direct measurement of the 
amplitudes of the corresponding power-laws). Thus no 
parameter fitting is applied. We observe an excellent data 
collapse for the entire range of the crossover confirming 
the phenomenological ansatz. In the inset of Fig.[21 we 
plot the same data without metric factors. As can be seen 
each model is characterized by its own scaling function. 

Since the entire crossover region covered several 
decades it could be difficult to observe small but sys- 
tematic differences between the scaling functions of both 
models. It its therefore instructive to examine the 
crossover via the so-called effective exponent 
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(12) 



The corresponding data are shown in Fig.|21 The ex- 
cellent data collapse of of both models over more 
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FIG. 2; The rescaled order parameter. The metric fac- 
tors are given by c = apO^ and d = Or. The data of 
both models display an excellent collapse to the universal 
crossover scaling function ^}\{x, 1). The dashed lines cor- 
respond to the asymptotic behavior of the two-dimensional 
system {/3n=2 = 0.639 HI) and of the mean-field behavior 
(/3mf ~ 1). Neglecting the metric factors each model is char- 
acterized by its own scaling function (see lower left inset). 
The upper right inset displays the effective exponents /3cft for 
both considered models. The data of the CLG model and of 
the CTTP exhibit an excellent data collapse. 



than 7 decades strongly supports the hypothesis that the 
crossover function is a universal function. 

We now consider the order parameter fluctuations. 
Analogous to the order parameter we make the scaling 
ansatz (7;^^ = 0) 

a^Ap^{p,R^,,) ^ S(a,(p-p,,„) A,a^ii?-/ A^). (13) 

Again the mean-field behavior should be recovered for 
R ^ 00, implying t){x,0) — Dtap{x,0) — 2 as well as 
Qa = aA,H^oo- Setting a~^R~J X'^ — 1 yields 

a^Ap^{p,R^,,) ^ S)(a,(p-p.,«) a^,/^i^i^^l). 



(14) 



For finite R, the fluctuations diverge at the critical point, 
i.e., the universal function J) scales as 



2)(a;,l) 



for 



0. 



(15) 



The universal amplitude to^ ^ can be determined in the 
following way: The scaling form Eq. (|13|l has to equal for 
R = 1 the D-dimensional scaling behavior [see Eq. 



Ap. ~ 0-1^^,^(1,0) 



jP - P..R=l) 
Pc,R=l 



-7d 



(16) 



Thus we find 



TOa.p = ^(1,0)' 



Pc,i?.= l ^p,R.- 



^p,R=l Pc,R- 



(17) 
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FIG. 3: The rescaled fluctuations of the order parameter. 
The metric factor is given by c = OpOR. The data of both 
models display s good collapse to the universal crossover scal- 
ing function T){x,l). The dashed lines correspond to the 
asymptotic behavior of the two-dimensional system (7^=2 = 
0.381 13] and rriA.p = 1.28) and of the mean-fleld behavior 
(Tmf = 0). The inset displays the corresponding effective 
exponent 7^ff. 



where the value of the universal scaling function 
D{1, 0) = 1.87±0.11 is obtained via direct measurements 
of the corresponding two-dimensional systems. Accord- 
ing to the scaling form Eq. I|14() we plot in Fig.O the 
rescaled fluctuations as a function of the rescaled control 
parameter for the two-dimensional CLG model as well as 
for the CTTP. We observe again a good collapse of the 
data over the entire region of the crossover. Furthermore, 
both asymptotic behaviors are recovered, confirming the 
scaling ansatz Eq. H13() . 

The corresponding effective exponent 7'jj. = 
d\n'I){x,l)/d\nx is displayed in the inset of Fig.|31 
Although the data of the effective exponent are suffering 
from statistical fluctuations one can see that both models 
are characterized by the same universal behavior. 

At the end we consider the critical amplitudes of the 
scaling functions. Using the above discussed scaling 
forms it easy to show that these amplitudes display a 
singular dependence on the range of interactions. For in- 
stance the order parameter scales sufficiently close to the 
transition point {x 0) as [see Eqs. (0 110lll|l ] 



The inset of Fig.^ shows the corresponding data for the 
CTTP. As can be seen, the data of various interaction 
ranges tend to the same power-law behavior if one ap- 
proaches the transition point. 

In conclusion, the crossover from mean-field to non- 
mean-field scaling behavior is numerically investigated 
for two different models exhibiting an absorbing phase 
transition. Increasing the range of interactions we are 
able to cover the full crossover region which spans sev- 
eral decades of the control parameter. The excellent col- 
lapse of the effective exponents of both models strongly 
supports the interpretation of the crossover scaling func- 
tions in terms of universality, i.e., the crossover function 
is universal. 
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ni0D-f3MF)/<t> 



Pc,R=l 



(18) 



Thus this scaling law and the corresponding scaling law 
for the fluctuations are only valid for finite interaction 
ranges whereas they become useless for infinite R, sig- 
naling the change in the universality class for i? — > cxd. 
This amplitude scaling can be observed in simulations. 



TABLE II: The non-universal metric factors determined from 
previous simulations via direct measurements of the corre- 
sponding power-laws. 



pc,R=l ap^R = i flA,i? = l Pc,R- 



Cip,R^o 



CLGi,=2 0.34494 0.5089 15.50 0.1244 0.1635 12.02 
CTTPd=2 0.69392 0.3410 50.18 1/2 0.3345 24.85 



